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ABSTRACT 

We address the question of whether the distribution of warm-dust compositions in IR- 
bright galactic disk PNe (paper I, Casassus et al. 2000), can be hnked to the underlying 
stellar population. The PNe with warm dust emission represent a homogeneous popu- 
lation, which is presumably young and minimally affected by a possible dependence of 
PN lifetime on progenitor mass. The sample in paper I thus allows testing the predic- 
tions of single star evolution, through a comparison with synthetic distributions and 
under the assumption that tip-of-the-AGB and PN statistics are similar. We construct 
a schematic model for AGB evolution (adapted from Groenewegen & de Jong 1993), 
whose free-parameters are calibrated with the luminosity function (LF) of C stars in 
the LMC, the initial-final mass relation, and the range of PN compositions. The ob- 
served metallicity gradient and distribution of star forming regions with galactocentric 
radius (Bronfman et al. 2000) allow us to synthesise the galactic disk PN progenitor 
population. We find the fraction of 0-rich PNe, f(0), is a tight constraint on AGB 
parameters. For our best model, a minimum PN progenitor mass M™™ = 1 Mq pre- 
dicts that about 50% of all young PNe should be 0-rich, compared to an observed 
fraction of 22%; thus M™" = 1.2 M©, at a 2a confidence level (M™" = 1.3 M© at 
la). By contrast, current AGB models for single stars can account neither for the 
continuous range of N enrichment (Leisy & Dennefeld 1996), nor for the observation 
that the majority of very C-rich PNe have Peimbert type I (paper I). f(0) is thus an 
observable much easier to model. The decrease in f(0) with galactocentric radius, as 
reported in paper I, is a strong property of the synthetic distribution, independent of 
M™'". This trend reflects the sensitivity of the surface temperature of AGB stars and 
of the core mass at the first thermal pulse to the galactic metallicity gradient. 



Key words: planetary nebulae: general - infrared: ISM: lines and bands - ISM: 
abundances - stars: evolution - stars: AGB and post- AGB. 



1 INTRODUCTION 

The vertical stratification of planetary nebulae (PNe) 
according to their N/0 abundance ratio (i.e. the Peim- 
bert 1978 classification), is an indication that PN com- 
positions can be linked to the progenitors' masses. Yet 
few attempts have been taken at building a comprehen- 
sive picture of the distribution of PN chemical compo- 
sitions in the galactic context, although the range of 
C, N and O abundances has been successfully modelled 
(e.g. Kaler et al. 1978, Renzini & Voh 1981, Bryan et 
al. 1990 , Groenewegen et al. 1995, Henry et al. 2000). 
The difficulties in PN population synthesis stem from 
an incomplete understanding of the late stages of stellar 
evolution, even for single stars on the asymptotic giant 



branch (AGB). Additionally, knowledge of the observa- 
tional properties of PNe suffer from the lack of accurate 
distances, and the variety of object types and ages. 

In paper I (Casassus et al. 2000) we reported on 
the 8-13 /^m dust features in catalogued compact and 
IR-bright galactic disk PNe. The PNe with warm dust 
emission represent a homogeneous population: their 
maximum age is limited by the presence of dust close to 
the central star. The sample is thus minimally affected 
by a possible dependence of PN lifetime on progenitor 
mass (keeping in mind that selection effects in the direc- 
tion of differential 'warm-dust' lifetimes as a function of 
grain composition are difficult to quantify, see Section 6 
in paper I and Section ^ in this work) . The dust-grain 
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spectral signatures also provide a diagnostic of the C/0 
chemical balance, and an alternative to the gas-phase 
abundances (e.g. Barlow 1983, Roche 1989). There is a 
trend for a decreasing fraction of 0-rich PNe outside 
the solar circle, reflecting the decrease in the ratio of M 
to C stars found by Thronson ct al. (1987) and Jura et 
al. (1989). 

Can the current scenarios for single star evolu- 
tion reproduce the observed frequencies of O- or C-rich 
grains in PNe, and the frequencies of Peimbert (1978) 
types among grain-types? Does the metallicity depen- 
dence of AGB evolution models predict the trends with 
galactocentric radius? In this article, we analyse the 
galactic distribution of PNe compositions in terms of 
simple stellar evolution models. We employ a set of ana- 
lytical prescriptions in order to evolve stars with a con- 
tinuous range of initial masses (Afi, from 1 to 7 Mq) 
and initial metallicity {Zi, from 0.005 to 0.035). This 
synthetic AGB model, coupled with a simple model for 
the galactic disk, allows us to check whether single star 
evolution can reproduce the distribution of PN chemical 
abundances. 

We will conclude that the fractions of PNe display- 
ing the dust signatures of 0-rich and C-rich grains, 
f(0) -I- f(C) = 1, are tight constraints on the param- 
eters of AGB evolution models. Synthetic AGB models 
that account for f(0) also reproduce the C star lumi- 
nosity function (LP) in the LMC; but the reverse is not 
true. However, for the PN sample analysed here, the 
minimum progenitor mass needs to be adjusted - our 
analysis favours a minimum progenitor mass of 1.2 M© 
at 2a. The galactic f(0) gradient stems from the metal- 
licity gradient, through the prescriptions for the core 
mass at the first thermal pulse and the AGB surface 
temperature. These results are an incentive to perfect 
the modelling, and pursue the use of PNe as probes of 
their progenitor population in different galactic envi- 
ronments. Also, the advantage of PNe for testing AGB 
evolution models resides in their well defined evolution- 
ary status. For comparison, it is difficult to ascertain the 
TP-AGB membership of M stars, making the trends in 
the C/M star ratio much more difficult to model (see 



Section [A2T| ). 

In Section we present the procedure followed to 
synthesise the galactic PN population, together with 
a brief description of the model used to evolve stars 
along the AGB (complete details are given in the Ap- 
pendix). The free parameters in the AGB model were 
calibrated with observed quantities; this synthetic AGB 
model should thus be taken as an observational require- 
ment rather than a detailed calculation. We investigate 
whether a model that reproduces the C star LF in the 
LMC also reproduces f(0). Section y describes the pro- 
cess of linking the model results with the PN sample, 
the results of which will be presented in Sect. ^ to- 
gether with an estimate of the minimum mass required 
to model the population of PNe with warm dust emis- 



sion features. Sect, a contains a summary of our con- 
clusions. 



2 SYNTHETIC DISTRIBUTION OF PNE 
COMPOSITIONS IN THE GALACTIC 
DISK 

2.1 Synthetic AGB model 

We aim at testing existing AGB evolution scenarios, 
such as that proposed by Groenewegen & de Jong (1993, 
GJ93). Synthetic AGB evolution models use empirical 
laws derived from detailed stellar evolution codes. As a 
starting point, we used GJ93, who follow Renzini and 
Voli (1981, RV81) and include an explicit metallicity 
dependence whenever possible. Other models exist, like 
those of Bryan et al. (1990) and Marigo et al. (1996, 
1999), but they have not been tested as extensively 
against observations as that of GJ93, which is also very 
light on computer resources and readily reproducible. 
The model is primarily concerned with the thermally 
pulsing AGB (TP-AGB) , while the initial conditions at 
the start of the TP-AGB are functions of {Mi, Zi). 

The set of analytical prescriptions was taken from 
the compilation by GJ93, although we also considered a 
different treatment of hot bottom burning (HBB), mak- 
ing use of the results in Forestini & Charbonel (1997, 
FC97)[|. The free parameters which bear directly on 
f(0) are 77agb, the mass loss parameter on the TP-AGB 
for a scaled Reimers (1975) law, and A, the fraction 
of core growth between two consecutive thermal pulse 
that is dredged-up and mixed. We kept fixed the other 
parameters in GJ93: the minimum core mass for 3'^^ 
dredge-up, M™™ = 0.58 Mq; the mass loss parameter 
for the red giant branch, t^rgb — 0.86; the mass loss 
parameter for the early- AGB, t^eagb = 3. 

The model was required to reproduce the C star 
LF in the LMC and the initial-final mass relation for 
stars in the solar neighbourhood, as in GJ93 (Sect.|A2[). 
The best model, which fits the C star LFs from Costa 
& Frogel (1996) and GJ93, corresponds to A = 0.75, 
'7AGB = 4, with a mass-loss rate high enough to keep 
the final masses within the observed limits. However, 
while an equally good fit to the GJ93 C star LF is ob- 
tained with A = 0.6, ?7agb = 5, this is rejected by f(0) 



(Section [4.21). 

The resulting PN abundances (A = 0.75, 77agb ~ 
4) , as a function of initial mass and initial metallicity are 
shown in Figure |l|. The set of abundances summarised 
in Fig. 1^ will be referred to as our standard model. The 
stellar ejecta were averaged over the last 25000 yr, to 
compare with gas-phase abundances. The different PN 



* In GJ93, the density at the bottom of the convective envelope 
was taken as a constant 2 g cni~^ (Martin Groenewegen, private 
communication), while these values are reached only in the most 
massive models of FC97 
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classes based on the relative abundances of ^^C, ^^C, 
N, O have their niche: O rich PNe are expected for 
both the lowest and highest Mi (type I PNe occurring 
for the highest progenitor masses), and ^^C/^'^C ranges 
from CN equilibrium values (^^3) to >100. 

If the stellar ejecta is averaged over the last 2000 yr 
only, which is taken to represent warm-dust compo- 
sitions (see Section 0), the contours for PN compo- 
sitions are still the same as in Fig. H in the case of 
log(N/0) (albeit a noisier Ila contour). In the cases 
of C/0 and ^^C/^'^C the PN composition map as a 
function of {Mi,Zi) is also similar to Fig. |l], but the 
point to point variations are sharper and reach higher 
peak value ('^13 and 3000 at tpN = 2000 yr against 
~7 and 900 at ipN = 25000 yr, for C/0 and i^C/i^C, 
respectively). There is one important change between 
ipN = 2000 yr and ipN = 25000 yr, however: in terms 
of the PN compositions, for tp^ = 2000 yr the effects 
of HBB for masses higher than 4-5 Mq are not as dra- 
matic, and C/0> 1 even for the most massive progen- 
itors (see Frost et al. 1998 for a physical discussion on 
the quenching of HBB). 

2.2 Problems with N enrichment in PNe 

In order to reproduce the range of N and C enrich- 
ment found in PNe and to account for the observed 
frequencies of type I and Ila PNe, we find that HBB 
must be much more efficient than currently modelled. 
As noted by Leisy & Dennefeld (1996), there is conti- 
nuity in the nitrogen abundances of PNe in the LMC, 
with no particular grouping above the type I thresh- 
old. And yet at Z=0.008, typical of the LMC, the syn- 
thetic AGB model predicts no type Ha PNe at all (in 
the Faiindez-Abans & Maciel 1987 classification, i.e. 
-0.6 < log(N/0) < -0.3). These conclusions also hold 
for the GJ93 model: with their treatment of HBB, Fig. [^ 
is largely unchanged (except perhaps at Zi ~ 0.03, 
where the type I threshold is at 4.6 Mq instead of 
5.5 Mq). Models in which core-growth on the AGB is 
negligible, and where HBB effectively occurs above a 
critical core mass, predict a sharp transition in N/0 
ratio and progenitor mass, which is not observed. 

Additionally, the majority of very C-rich PNe 
(those whose mid-IR spectra display the family of emis- 
sion bands associated with PAHs) are of type I (pa- 
per I). However, the existing synthetic AGB models can- 
not produce N-rich and very C-rich PNe in the observed 
frequency. Fig. n^ shows how the highest C enrichment 
should be found for PNe with solar N/0 ratios (again, 
the PN composition map obtained with a 2000 yr aver- 
age of AGB mass- loss leaves Fig. |l| largely unchanged). 
In terms of the gas-phase abundances, type I PNe are 
divided in their C/0 chemical balance (Leisy & Den- 
nefeld 1996, their Fig. 12; see also AUer & Czyzak 1983 
for galactic disk PNe) . But the sharp onset of HBB with 
Mi produces 0-rich PNe, exclusively (for Zi > 0.01). 
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Figure 1: PN abundances averaged over the last 
25 000 yr of AGB evolution as a function of initial 
mass in x, initial metallicity in y. From top to bottom: 
log(N/0), with darker grey scale contours at -0.8, - 
0.7, ..., 0.6, and the thin solid lines at -0.6 and -0.3 
represent the contours for type Ha and type I PNe, 
while the thick line traces the minimum initial mass 
for the occurrence of the second dredge-up (a sharp 
transition at —0.3 results in missing contour levels, but 
the minimum and maximum plotted values of the grey 
scale are -0.8 and 0.6); C/O, with darker contours at 
log(C/0)= -0.6, ..., 0.4, 0.6, and the solid line is at 
C/0=1; 12C/13C, darker contours at log(i2C/i3C)= 
0.7, 0.9, ..., 2.9, and the solid line is at i2c/i3c=30. 
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An exception to the above problems with synthetic 
AGB models is provided by RV81. It appears that RV81 
underestimated the mass loss rate (??agb ~ 1/3 predicts 
excessive final masses, GJ93 and references therein), 
and the extended lifetime allows the temperature at the 
bottom of the convective envelope to rise significantly 
with the growth of the core during the TP-AGB. A 
3.3 M0 star can start dredging-up carbon and then un- 
dergo HBB later in its evolution, in a balance that pro- 
duces C and N rich PNe. Thus the problems with C and 
N enrichment in PNe could be solved by artificially fix- 
ing the slope for the increase of HBB temperature with 
core-growth, dTB/dMc, to the values found by FC97 
for their 5-6 M© models. We found that models with 
dTB/dMc = 2 10^ K Mq^ produce type I and Ila PNc in 
the observed frequencies (around 30% and 20% in the 
sample of paper I, see also Maciel &: Dutra 1992), with 
type I PNe being preferentially very C-rich. An interest- 
ing alternative to a high dTs/dMc slope can be found 
in Siess & Livio (1999). The accretion of brown dwarfs 
or giant planets could trigger CNO processing in the 
convective envelope. Although the concept of envelope 
heating from the kinetic energy of accreted material is 
appealing, it is not obviously related to the observed 
N/0 latitude effect. 

Thus the N/0 classification of PNe, perhaps one of 
their most widely used chemical characteristic, turned 
out to be impossible to use in PN population synthesis, 
at least with the models we tried. On the other hand, 
f(0) depends mostly on M™'" coupled with the core 
mass at the first thermal pulse (in this formalism). In 
fact, above a threshold initial mass, all PNe are pre- 
dicted to generate C-rich warm dust during the earliest 
stages of their evolution. f(0) is therefore a much easier 
quantity to model. 



2.3 Synthetic distribution of PN progenitors 
in the galactic disk 

A set of 10000 stars was evolved to construct the syn- 
thetic distribution of PNe in the galactic disk. The 
probability distribution for the star formation rate as a 
fimction of galactocentric radius was extrapolated from 
the ultra-compact H 11 region survey in Bronfman et al. 
(2000). We assumed the number density of massive star 
forming regions would also trace the number density of 
newly born PN progenitors. Initial stellar masses were 
distributed according to an initial mass function (IMF) 
index of —1.72 (in dN/dlogM). This index is also used 



in the derivation of the LMC C star LF (Section A2.1), 



and is the same as in GJ93. The lower mass cutoff for 
the compact PNe progenitors was kept as a free param- 
eter, M™™, and the upper mass limit was 7Mq. The 
IMF index may seem steep for intermediate mass stars, 
and we investigated using —0.95, as derived in Sabas 
(1997) from a complete sample of B5-F5 stars and hip- 



parcos parallaxes (over 1.2-4 Mq). In this case M™™ 
turned out to be somewhat lower (see Section ^). 

Radial diffusion in the galactic disk was taken into 
account following the results presented in Koppen and 
Cuisinier (1997), in the epicyclic approximation. The 
age- velocity relation by Wielen (1977), 



/l00 + 600^^kms- 
10^ yr 



(1) 



in conjunction with a velocity ellipsoid of constant 
shape (Wielen 1977), 



cTu : c^v : o-w = 0.79 : 0.46 : 0.41, 



(2) 



allows to estimate the spread of the normal distribu- 
tion by which final galactocentric radii are distributed 
(Wielen et al. 1996), 

crR(age) = 1.56cru(age)/K, (3) 

where k = 31.6kms~^ kpc~^ is the epicyclic frequency 
of stellar orbits in the disk, at R = Rq. The factor 1.56 
could altogether be neglected in this schematic treat- 
ment, which will be applied to ~ .5 — 1.5 -Ro- 

As argued in Carraro et al. (1997), we adopted a 
radial metallicity gradient of —0.07 dex/kpc and con- 
stant in time, and assumed a 0.005 FWHM scatter when 
drawing initial metallicities. Z(R = Ro,t) was required 
to match the age metallicity relation by Meusinger et 
al. (1991), 



[Fe/H](i) « log[Z{R = Ro,t)/ZQ] = 

log [1.98- 1.88(1 -t/28.5)i-2^] 



(4) 



3 LINKING THE MODELLED AND 
OBSERVED PN GALACTIC 
DISTRIBUTIONS 

We based the comparison between modelled and ob- 
served distributions on the sample of compact and IR 
bright PNe from paper I. In linking the modelled and 
observed PN galactic distributions, the selection effects 
in the PN catalogues need to be taken into account, 
as well as the exact relation between the AGB stellar 
ejecta and the observed PN compositions. 

The PNe in the sample with warm dust emission 
were required to have IRAS 12 /im band fluxes in ex- 
cess of 0.5 Jy. Their distribution was modelled by ap- 
plying the same selection criterion to the synthetic PN 
population. It is shown in paper I that the PNe with 
warm dust should be optically thick in the Lyman con- 
tinuum, in which case the four IRAS bands give a good 
representation of their bolometric flux, with a fraction 
of 0.25±0.14 of the total fiux in the 12^m band. The IR- 
bright selection criterion should probe the whole of the 
galactic disk. But selection effects that limit the galactic 
disk sampling area come from the optical discovery sur- 
veys. Extinction in the galactic disk needs to be taken 
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into account. Condon & Kaplan (1999) report extinc- 
tion values towards 429 PNe; no PNe are found with 
c — log{iif3'^°'^^ /lif3°^^) in excess of 3.3, with an approxi- 
mate cutoff in the distribution at c ~ 2. They also show 
that c depends strongly on galactic coordinates, with 
values of order 1 towards the galactic anticentre. The 
bulk of the extinction to PNe is thus galactic, rather 
than intrinsic, and a maximum value of c=2~3.5 would 
give ^v = 3.8 — 6.6. We adopted a maximum interstel- 
lar extinction of Ay = 4, as for the objects with warm 
dust intrinsic extinctions of at least one magnitude have 
been reported. 

A gaussian profile|j in galactocentric radius for 
A(i?), the 'linear density of extinction', gives 



AyiDJ) 



I'D I'D 

/ dsA{R) = / dsA exp(- 
Jq Jo 



2Rl 



)• (5) 



The values oi A = 4 magkpc^^, Ra = 3.6 kpc were 
fixed by requiring Aqc = 20 mag for the extinction 
to the galactic centre, and a minimum galactocen- 
tric radius of 0.5 -Ro for Ay'^^—A. These requirements 
give a local linear density of visual extinction of only 
0.25 magkpc"^, and the value of Aqc — 20 mag may 
seem low as well. However, we did not take into account 
the dependence of A{R) on height above the galactic 
plane, and that PNe would preferentially be discov- 
ered towards directions of low extinction: if the aver- 
age linear density of extinction to local PNe were really 
1 magkpc"^, no PNe would be detected beyond 4 kpc. 
It is apparent from Fig. 6 in paper I that the maximum 
distance to galactic disk PNe with warm dust is about 
8 kpc, which we adopted as a maximum distance in the 
synthetic population. This estimate of the sampled area 
in the galactic disk is intented to approximate the ob- 
served face-on distribution. 

In the case of the PNe with warm dust emission 
features, the problem of a lifetime dependence on the 
progenitor mass is greatly simplified by the presumed 
youth of the objects. We assumed that the compact and 
IR-bright PNe with warm dust emission all have the 
same lifetime, neglecting any dependence on Mi. But it 
should be kept in mind that the durations of the 'warm- 
dust' phase lifetimes could depend on the dust opacities. 
It was suggested to us (Kevin Volk, private communica- 
tion) that this could be the reason why the fraction of 
0-rich objects is f(0)~22% from the dust grain compo- 
sitions, while it is closer to 40% from the plasma diag- 
nostics. Some arguments against this possibility may be 
found in Section 6 of paper I, and the next section de- 
velops the idea that the difference in f(0) could be due 



T Koppen & Vergely (1998) modelled the galactic disk extinction 
by matching an exponential linear density of extinction in R and 
z to bulge PNe. But the dust and molecular content at i? < 2 is 
similar to that of the molecular ring at i? '^ 0.5 Rq (Bronfman et 
al. 1988, Deul & Burton 1991), and it is not clear whether there 
is a sharp central concentration of dust in the galactic disk. 



to the different timescale over which the AGB ejecta 
should be averaged (ipN in this work). In any case, the 
conclusions we reach on f(0) being a tight constraint 
on AGB parameters are valid even if f(O)~50%. 

A lower metallicity limit for the production of dust 
is discarded, as even some very low metallicity PNe 
show warm-dust emission (e.g. in a forthcoming article 
on bulge PNe, we will report on the type IV PN M 2-29, 
that shows the silicate signature in spite of a very low O 
abundance of [0/H]= —1.4 dex, as estimated by Ratag 
et al. 1997). 

Since PNe with no warm dust correspond to later 
evolutionary stages (paper I), the link between the com- 
puted C/0 ratios and the dust features is very tight. As 
discussed in paper I, Silicates represent C/0<1, and 
SiC and UIR bands C/0>1. The only problem left 
in linking the synthetic and observed PN populations 
is defining how to compute the average compositions. 
If the warm dust grains are heated directly by radia- 
tion from the central stars, then a temperature of 100- 
200 K places the grains at ~ 10^^ cm from the exciting 
star, which corresponds to an ejection age of 400 yr at 
lOkms^^. In this picture, the warm dust grain composi- 
tion derives from the very last stages of AGB evolution. 
Thus the PNe compositions were averaged over the last 
2000 yr, to compare with the dust emission features, 
while a more topical 25 000 yr was used to compare 
with the gas phase abundances. 

The resulting synthetic distribution can be seen in 
Figure 0, in the case of a minimum progenitor mass 
of 1.2 M0. Galactocentric azimuths were drawn from a 
uniform random distribution. 



4 RESULTS: THE MINIMUM MASS FOR 
PN PROGENITORS 

4.1 The minimum mass for PN progenitors 

Having calibrated all the free-parameters in the syn- 
thetic AGB model, the mass range for PN progenitors 
still needs to be specified. The lower mass limit has 
a direct influence on the synthesised PN distribution, 
although the exact upper mass limit is of secondary im- 
portance, high mass stars being less frequent. 

The property of the synthetic distribution that can 
be most directly compared to the observed distribu- 
tion of PNe with warm dust is the relative frequency 
of C-rich and 0-rich nebulae. Shown in Figure H is the 
fraction of O rich PNe as a function of minimum pro- 
genitor mass, Af™'", for the set of parameters given in 
PN compositions were computed by aver- 



Section 2.1 



aging the stellar mass loss over the last tpN = 2000 yr 
of evolution on the AGB to compare with the dust sig- 
natures. The shaded areas correspond to the range of 
observed values within ± 1 cr (from paper I) , 



f(Orich) = 0.22 ±0.06. 



(6) 
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Figure 2: Synthetic distribution of PNe in the galac- 
tic disk, with filled circles representing 0-rich nebulae, 
while open circles represent C-rich nebulae. The dotted 
lines correspond to the solar circle (Ro=8.5 kpc) and to 
1 = 0. The clear area is the region where sources would 
have been included in the PN catalogues, within the 
adopted limit of Ay = 4 and dmax = 8 kpc. The AGB 
mass loss was averaged over the last 2000 yr of evolu- 
tion. For the clarity of the plot, the initial number of 
PNe is only 1000, and comes down to ^200 within the 
galactic disk sampling area. This Figure can be com- 
pared to Fig. 6 in paper I. 
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Figure 3: The fraction of 0-rich PNe, f(0), as a function 
of minimum main sequence progenitor mass, M™'". The 
shaded area correspond to the observed values within 
la. 

In order to match the observations, Af ™'" is constrained 
to 



Af" 



> 1.2M 







(7) 



at two-cr confidence level (or ^1.3 M0 at one-cr). The 
IMF slope in our standard model is rather steep, and 
for comparison a shallow IMF of index —0.95 (Sabas 
1997), gives Af""" > 1.2 M© at one-a - f(0) does not 
seem very sensitive on the IMF. 

We will now test the value of A/""" against less 
accurate properties of the observed PN distribution, 
and analyse trends with galactocentric radius. In the 
synthetic PN distribution, the relative number of PNe 
inside and outside the solar circle is Nin/Nout = 1-3 
(rather independent of Ai^™'"), while the observed ratio 
is 0.9-1.6 (from Tables 6 and 7 in paper I), depending 
on the PN distance scale. The approximate agreement 
is an indication that the galactic disk sampling area of 
Section is appropriately defined. 

Figure ^ shows the fraction of O rich PNe, f(0), 
for R < Ro and R > Rq. There is a clear decrease in 
the model f(0), for any value of A/™'", and it is thus a 
robust property of the synthetic distribution. From pa- 
per I, the observed f(0) is 30 ±9% inside the solar circle, 
and 14 ± 7% outside. Figure Qd shows the proportion of 
type I PNe as a function of minimum progenitor mass. 
For any A/'°'", the fraction of type I PNe is far too low. 
This result can also be reached by estimating the pro- 
portion of progenitors between Af ™'" and ~ 4 Mq , even 
for a shallow IMF. Either this is yet another deficiency 
of current models for N enrichment, or the PN sam- 
ple is dramatically affected by biases against low mass 
progenitors. The latter case seems rather contrived; a 
complete PN sample is a working hypotheses for this ar- 
ticle, and in the event of a strongly biased sample, Af •"'" 
represents a lower mass limit above which we trust the 
observed f(0) - as well as f(type I). Additionally, we 
know N enrichment is not satisfactorily modelled (e.g. 
Section^. 
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Figure 4: Properties of the synthetic PN population as 
a function of minimum main sequence progenitor mass, 
M™™. a) f(0) inside and outside the solar circle, b) the 
fraction of type I PNe (the observed values is 29 ± 6%, 
from paper I and references therein). 



For comparison the gas phase abundance distribu- 
tion synthesised with tpN — 25000 yr predicts 5-12% 
more 0-rich PNe, depending on M""", and around 
28% for M™'" = 1.4 Mq (which gives f(0)=22% if 
ipN — 2000). This is due to a contribution of 0-rich 
objects produced by HBB, which is quenched over the 
last 2000 yr. In the sample of paper I, this fraction is 
40 ± 8%, and is - 40% in Zuckerman and AUer (1986, 
note that their sample is inhomogeneous, with a broad 
range of PN ages). The gas phase C/0 ratio is available 
for an even more restricted number of sources (of the 
ones with warm dust), and are subject to large uncer- 
tainties stemming from the difficulties in measuring the 
C abundance, especially in compact and IR-bright PNe 
- compare for instance AUer & Czyzak (1983) and Henry 
et al. (2000), in the cases of NGC 3242 and NGC 6826. 
Considering such uncertainties, we believe the difference 
between the fraction of 0-rich PNe as inferred from the 
dust compositions and the plasma diagnostics is due to 

tpN- 

In the framework of an artificially high dT-Q/dMc 
(Section [2^), the eS'ects of fpN on f(0) are stronger be- 
cause HBB occurs for lower progenitor masses and is 
thus more frequent. Accounting for the observed pro- 
portions of type I and Ha PNe (29 ± 6% and 12 ± 6%, 
paper I) would give an upper limit to M™'" of about 
1.4 Mq (above this limit, there are too many type I 
PNe). With M™'" = 1.2 M© and ipN = 2000 yr, the 
fraction of C-rich PNe that have Peimbert type I is 35%, 
and the observed proportion of type I PNe among PNe 
with C-rich dust grains is 37± 10% (paper I). Similarly, 
in the synthetic distribution the fraction of 0-rich PNe 
that have Peimbert type I is 26%, and the observed pro- 
portion of type I PNe among those with 0-rich dust is 
22 ± 14%. 



4.2 Model sensitivity 

What is the model sensitivity of the predicted frac- 
tion of 0-rich PNe? Of the relations used in the syn- 
thetic AGB model, which are the ones responsible for 
the galactic gradient in f(0), i.e. which are the most 
sensitive on the galactic metallicity gradient? 

The predicted fraction of 0-rich PNe is a sensitive 
function of the free parameters, as shown in Fig. pK. 
The sets of parameters correspond to a range of 3'^'^ 
dredge-up efficiencies, over an AGB star's lifetime (an 
extended lifetime results in more dredge- up episodes). 
But, at first sight of Figs, ob and rac, all models seem 
to acceptably reproduce the C star LF used in GJ93 
or that from Costa & Frogel (1996, CF96, with la un- 
certainties of ~ 10% at the peak). In terms of the C 
star LF x^ values, which we calculated summing over 
each logarithmic luminosity bin, from 4 to 5.8 for the 
CF96 LF and from 4.125 to 5.625 for the GJ93 LF, 
the AGB parameters adopted by GJ93 (?7agb = 4 and 
A — 0.75) seem reasonable. We assumed the GJ93 LF 
was built with 298 sources, about three times less than 
in CF96 (895 sources). But the associated x^ confidence 
levels are very low (with 7 and 10 degrees of freedom 
for GJ93 and CF96, respectively). It is difficult to as- 
sess the quality of the C star LF fit, from one set of 
parameters to another. It turns out the C star LF x^ 
values are very noisy: for the set of AGB parameters 
(A,77agb) = {(1,1), (0.9,2), (0.9,3), (0.6,3), (0.75,4), 
(0.6,4), (0.75,5), (0.6,5), (0.5, 6)}, x' =(11, 6, 14, 7, 6, 
8, 10, 6, 11) for the GJ93 LF and x^ =(41, 39, 24, 56, 
32, 39, 28, 27, 40) for the CF96 LF. On the other hand, 
f(0) exhibits broad variations over the same set of AGB 
parameters, and the last two models with 77AGB — 5 are 
clearly rejected. Models with ryAGB — 2 and A < 0.6 



predict excessive final masses (see Section A2.2), and 
j^vam < j^ js unlikely considering the extremely low 
proportion of type I PNe (around 3%). We stress that 
the constraints put on the AGB parameters are model- 
dependent, rejected parameters could well fit another 
set of AGB prescriptions. Our purpose here is to illus- 
trate the usefulness of PN population synthesis, pro- 
vided a proper PN sample is built. 

Of the relations used in the synthetic AGB model 
(see Appendix), the following include an explicit metal- 
licity dependence: 

• the core-mass-interpulse period relation, 

• mass loss prior to the AGB, 

• core-mass at the first thermal pulse, 

• luminosity at the first thermal pulse, 

• core-mass-luminosity relation, 

• T,b{L,Z,M), 

• evolution rate, 

• hot bottom burning density and temperature. 

The synthetic PN population was tested for sensitiv- 
ity on the galactic metallicity gradient by fixing the 
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Figure 5: a) The fraction of 0-rich PNc in the solar neighbourhood, in sohd hues are model results from PN population 
synthesis as a function of minimum progenitor mass (in Mq). The coloured area corresponds to the observed fraction 
±1(7. r^AGB is the mass loss coefficient of a scalable Reimers law, and A is the 3''"^ dredge-up parameter, b) Synthesis 
(solid line) of the LMC C star LF used in GJ93 (dotted line) , for the same sets of parameters as in a), c) Synthesis 
(solid hue) of the observed (dotted hue) LMC C star LF from Costa & Frogel (1996). 
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metallicity entry of all the previously listed relations at 
Z = 0.02. The initial composition was kept as described 
in the Appendix, i.e. a metallicity-scaled solar mix, with 
a galactic metallicity gradient and an age-metallicity re- 
lation. The resulting synthetic PN population did not 
show a gradient in the fraction of 0-rich PNe. Subse- 
quently the constraint on the metallicity was relaxed 
one relation at a time, and the resulting model pop- 
ulation was inspected using the equivalent of Fig. |4|a. 
The difference in the fraction of 0-rich PNe inside and 
outside the solar circle was used to test the dependence 
upon the metallicity gradient. 

The most sensitive relations are Teff(L, Z, M) and 
the core-mass at the first thermal pulse. The other re- 
lations resulted in an insignificant change in the frac- 
tion of 0-rich PNe inside and outside Rq. The synthetic 
population generated using a flat metallicity profile was 
confirmed not to show a gradient in the fraction of O- 
rich PNe. The metallicity dependence of f(0) resides in 
the following: 

(i) Tcff increases with lower metallicity, which for a 
fixed luminosity translates into a lower envelope ra- 
dius, and a lower mass-loss rate. The extended AGB- 
lifetime thus leaves more time for core-growth and the 
C-enrichment of the envelope. The relation used here 
for Tcff, is a model result from Wood (1990, see also 
Vassihadis & Wood 1993). 

(ii) The core-mass at the first thermal pulse, Mc(l), 
is a decreasing function of metallicity. Together with 
a fixed minimum core-mass for the occurrence of the 
third dredge-up (the free parameter Af™'"), the lower 
limit progenitor ZAMS mass for the production of C- 
rich nebulae decreases with metallicity. 

A full search in all possible combinations in param- 
eter space is postponed to a future investigation, based 
on a more solid AGB model, and a larger PN sample. 



mum PN progenitor mass of 1 Mq predicts that about 
50% of all young PNe should be 0-rich, when only 22% 
is observed. The minimum PN progenitor mass required 
to reproduce the properties of the observed distribution 
of PNe with warm dust is Af™'" = 1.2, at a 2a confi- 
dence level. 

(ii) The observed decrease in the fraction of 0-rich 
PNe with galactocentric radius is a robust property of 
the synthetic distribution. 

(iii) Independently of Af '"'", values for the dredge up 
parameter A and the mass-loss rate parameter t^agb 
that reproduce the LMC C star LP, do not necessarily 
match the observed fraction of 0-rich PNe, f(0). 

The relative frequencies of C- and 0-rich PNe depend 
mostly on the value of the minimum core mass for C 
dredge-up, which is calibrated with the LMC C star LP 
(as in GJ93). But the treatment of HBB is of crucial 
importance in predicting PN compositions, as it directly 
accounts for the Peimbert (1978) types. Measurements 
of the ^^C to ^^C isotopic ratio in PNe can provide an 
observational test, since a transition from >100 to '^S 
is expected with the onset of HBB (e.g. RV81). 

Overall, the good agreement obtained with the ob- 
served PN distribution is an indication that the statis- 
tics of dust signatures in PNe can be used as probes of 
their progenitor population in various galactic environ- 
ments. 
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5 CONCLUSIONS 

PN population synthesis was investigated using a 
schematic model for AGB evolution, whose free- 
parameters were calibrated with the luminosity function 
of C stars in the LMC and the initial-final mass rela- 
tion. The galactic distribution of PN progenitors was 
generated using the observed metallicity gradient and 
distribution of star forming regions with galactocentric 
radius. The PNe with warm dust emission from paper I 
represent a homogeneous population, which is presum- 
ably young and thus minimally affected by a possible 
dependence of PN lifetime on progenitor mass. We ex- 
amined the consequences of assuming the statistics of 
PN compositions in the sample of paper I reflect tip-of- 
the-AGB values, and reached the following conclusions: 

(i) The fraction of PNe with 0-rich dust places con- 
straints on the range of PN progenitor masses. A mini- 



APPENDIX A: ANALYTICAL 
PRESCRIPTIONS FOR AGB EVOLUTION 
AND CALIBRATION OF FREE 
PARAMETERS 

Al Synthetic AGB evolution code 

To summarise, we will describe the steps followed in 
the code, for a star with fixed {Mi,Zi). All the ana- 
lytical prescriptions are written in GJ93, except when 
explicitly indicated as not being so. The initial com- 
position is calculated assuming the relative abundances 
of the elements are the same as in the solar mix. We 
follow the steps in GJ93 up to the first thermal pulse, 
taking the mass-loss prescription on the E-AGB from 
Groenewegen et al. (1995, their Eq. 1). The core mass 
and luminosity at the first thermal pulse, Afc(l) and 
L(l), are the initial conditions for the TP-AGB. Afc(l) 
is of significant importance, since it fixes for the most 
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part the initial-final mass relation, one of the obser- 
vational constraints. GJ93 take Mc(l) from Lattanzio 
(1989) for Mi under 3 Mq, and from RV81 above. The 
subsequent TP-AGB evolution is described by a loop, in 
which the interpulse period iip is determined by a core- 
mass-interpulse relation (from Boothroyd & Sackmann 
1988, GJ93 include a 'turn-on' phase). The luminosity 
is given by a core-mass-luminosity relation, with the fits 
in GJ93 (Boothroyd & Sackmann 1988 and Iben & Tru- 
ran 1978), who use a rectangular profile to account for 
the luminosity dip and flash. GJ93 calculate the rate of 
core-growth from the luminosity due to hydrogen burn- 
ing^ The mass loss rate is estimated via a tjagb scalable 
Reimers (1975) law. The stellar radius R is estimated 
from the (model) relation given by Wood (1990) be- 
tween the luminosity and the effective temperature of 
0-rich Mir as. The condition for ending the AGB is then 
tested (based on a critical envelope mass, Iben 1985), 
as well as the end of the interpulse period. In the case 
of a new thermal pulse, the condition for dredge-up is 
tested: Mc > M™'". M™'" is one of the free parameters 
in GJ93. 

The process following a third dredge-up event is 
the most sensitive part of this code, and it affects di- 
rectly the PN compositions. It is standard to model the 
dredged-up material enriched in He and ^^C as a frac- 
tion A '^ 0.6 of the core growth during the interpulse 
period AA/c, 



AMdredgod = AAMc. 



(Al) 



A could well be a function of metallicity or other pa- 
rameters, but is here taken as a constant. The com- 
position of the dredged-up material, the result of in- 
complete helium burning, is taken from Boothroyd & 
Sackmann (1988), Y^'^i-d" = 0.76, isc^'d-du ^ o.22, 
i6o3rd-du ^ 0.02. In the absence of HBB, the dredged- 
up material is simply diluted in the envelope. 

HBB is the main source of uncertainties in the 
model: existing AGB models for single stars have dif- 
ficulties in producing N and C rich PNe (see Section 
2.2). HBB was treated by setting the initial tempera- 
ture Tb(1) and density pb(1) at the bottom of the con- 
vective envelope to the values listed in FC97. We first 
interpolated in logTB(l) and logpB(l) as functions of 
Mc(l), for the two metallicity cases listed (LMC and 
solar), and then in logrB(l) and logpB(l) as functions 
of metallicity Z. For the subsequent evolution on the 
TP-AGB, we used the model results in FC97 as a func- 
tion of time to estimate dT-QJdMc and dpB/dMc- They 
obtained constant growth rates for pB and Tb, and we 
assumed core-growth to be roughly linear in time (as 
it is in GJ93), and given by Table 5 in FC97. A more 



accurate treatment is pointless, as major adjustments 
to dT^/dMc are required to explain the observations. 
Thus the results in FC97 led us to the following pre- 
scriptions for dp^/dMc and dT^/dMc, with upper limits 
dpB/dM^ < 10^ g cm-3 Mq\ Tb < 10^ K, pB < 1: 



log 



log 



dTB 
dMa 
dpB 

dMr 



8.33Mc(l) + 2.05, 
6.85A4(1) + 3.51, 


Z = 
Z = 


- 0.02, 
= 0.008, 


14.3Afc(l) - 10.4, 
13.1Afc(l) - 9.23, 


Z^ 

z^ 


= 0.02, 
= 0.008. 



(A2) 



(A3) 



As explained in Section E^, a significant increase 
in Tb during the AGB is required to explain C and N 
rich PNe. Equation A2 does not provide such a strong 
increase, except for the highest masses. As many un- 
certainties affect the treatment of HBB anyway, we 
investigated treating dTB/dM^ as a free parameter, 
and satisfactory results are obtained for dTB/dM^ = 
210^KMq\ 

Once Tb and pb are determined, dredged up 
material is exposed to CNO processing over a time 
^HBB, and material in the envelope is processed dur- 
ing iHBBiip/Tconv, whcrc Tconv ~0.5yr is the convec- 
tive turn-over time-scale|| (as quoted in FC97). We fixed 
^HBB = lO^^yrs so as to obtain ^^C/^'^C ratios at CN 
equilibrium values for the most massive stars, while 
avoiding the onset of the full tri-cycle (which is not 
reported to occur in the envelope of AGB stars, and 
reduces the O abundances to produce C/0 ratios of or- 
der 10). In the HBB treatment of GJ93, the envelope 
cannot reach ^^C/^'^C at equilibrium values, since they 
took only a very small fraction (^ 10^"*) of the envelope 
mass to be CNO-processed. In G J93 dredge-up material 
is fully exposed, with ^^C/^^C~3, but it is C-deficient, 
and ^'^C is diluted in the envelope. It is worth noting 
that even with these differences for HBB, the equivalent 
of Fig. n^ for the model in GJ93 HBB is identical - ex- 
cept for C/^^C. The injection of triple-alpha enriched 
material in the envelope is given by 



^ydu ^ Y 



3rd — du 



AAf, 



dredged; 



AX 



du 



1 



Xi{t)dt AMdrcdged, 



(A4) 
(A5) 



^BB 

with initial conditions for the abundance of element i, 
Xi{0) given as above, or otherwise taken as zero. The 
change in the surface abundances is 



Xf 



AXf"-f ^/Q™X(i)di 



A/c: 



AAfdrodgod 



(A6) 



where t^^B 



tuBB tip/ Tconv, and the initial condition 
for the time average of CNO processing is taken as 
Xi{0) = X°'^. The nuclear reaction rates for the CNO 



+ We did not include a luminosity increment due to HBB. The 
observational constraint is the C star LF, and a HBB correction 
would only affect massive M giants (and their mass loss ral^ 
which go in the C-phase at the very end of the AGB, see Fig. 



•a te') 

?E1 



3 The value of Tconv is of secondary importance. The total ex- 
posure time is tuBB iip/T"conv, and any change in Tconv can be 
compensated with tHBB ■ 
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tri-cycle were taken from Caughlan and Fowler (1988), 
and we followed the analytical treatment in Clayton 
(1968), with approximations for the branching ratios 
between the main CN-cycle and the two secondary cy- 
cles from Podsiadlowski (1989). 

As explained in Frost et al. (1998), there is a mini- 
mum envelope mass below which HBB is quenched. The 
critical envelope mass M^gg was adapted from GJ93, 



^W HBB - 


= 0.85Afp„, 


(A7) 


Mpn = 


= 6 (1.687- 9.092 Mc + 11.687M2 
-4.343 Af 3), 


(A8) 



where Mp^ is the RV81 criterion for PN ejection with 
b=0.5. The exact value of b and the 0.85 factor are of 
secondary importance, GJ93 having fixed these param- 
eters to match the results in RV81. 

We stress here that dredge-up must occur after 
HBB is quenched. Observational support can be found 
in Trams et al. (1999) and van Loon et al. (1999). It is 
of crucial importance for explaining the vast majority 
of N rich PNe with C-rich dust (corresponding approx- 
imately to the last 2000 yr of AGB mass loss) . 

A2 Observational constraints 

A2.1 The luminosity function of C stars in the LMC 

The population of stars in the LMC with initial masses 
between M and M + dM which are currently on the 
TP-AGB is 

N{M)dM (X SFR(M) x IMF(M) x tAGB{M)dM, (A9) 

where SFR(M) is the star formation rate in the LMC 
at time t(M), the time of birth as a function of mass, 
IMF(M) is the initial mass function, which will be 
taken as a power-law, and tAGB{M) is the lifetime 
on the TP-AGB as a function of main-sequence mass. 
There is an implicit metallicity dependence via the LMC 
age-metallicity relation. In practice, we started with 
N{M)dM ex SFR(M) x IMF(M), and weighted the re- 
sult with iAGB(-^^)- The procedure is equivalent to mul- 
tiplying each mass entry by tAGB,{M) (times an integer 
big enough to keep a satisfactory resolution) . 

The parameters describing the global evolution of 
the LMC were reproduced from GJ93. The star forma- 
tion rate in the LMC is approximated as SFR(M) ex 
exp(-i(M)/7), with t{M) in Gyr and an LMC age of 
TLMC=llGyr, the initial mass function is IMF(M) ex 
M~^-^^, and the age-metallicity relation is assumed to 
be linear in time, Z{t) = 0.01 - 0.008 x (1 - 1/13). The 
total lifetime as a function of mass, t{M), was taken 
from Iben & Laughlin (1989). 

Observational data with which to constrain the 
AGB population in the LMC is compiled in GJ93. The 
luminosity function of C stars is a very good constraint 
since there can be no doubt on their evolutionary status 
(i.e., at the luminosity implied by the LMC distance, C 




Figure Al: Time spent on the AGB (10^ years, in 
ordinates) as a function of initial mass, for the LMC 
model. The grey scale grows darker from the M to S to 
C phases. 

stars can only be on the TP-AGB). On the other hand 
the relative frequency of C and M stars is much more un- 
certain: the contamination from M stars on the E-AGB 
is difhcult to estimate. Hughes (1989) and Hughes & 
Wood (1990) found that most LPVs in the LMC have 
spectral type M5-f , and C/M= 0.63, provided all M 
stars on the TP-AGB are also LPVs. Based on previous 
measurements (Blanco & McCarthy 1983), GJ93 favour 
C/M> 0.63. But M giants spectral types are difficult to 
assess, and we feel it is safer to trust the C star counts, 
and admit that the number of M stars on the TP-AGB 
is not known, although probably greater than the num- 
ber of 0-rich LPVs - thus C/M< 0.6. With tu, is, ic, 
the time spent by a star in the M, S and C phases, the 
C/M giants ratio on the TP-AGB is given by 

JdMSFR{M) X IMF(M) x tc{M) 



C/M: 



(AID) 

Mo (we 



JdMSFR{M) X IMF(Af) x tuiM) 

where the integral runs from 0.93 M© to '^ 
kept the lower mass limit in GJ93). The M phase corre- 
sponds to C/O<0.81, the S phase to 0.81<C/O<l, and 
the C phase to C/0>1 (Smith & Lambert 1986). We fol- 
lowed GJ93 in neglecting obscuration of bright C stars 
due to excessive mass loss. GJ93 also estimated that not 
more than 3% of C stars brighter than Mboi = —6 could 
have been missed by the optical surveys. The functions 



^M, tg, tc are shown in Fig. Al. The formula above, as 
well as the Monte Carlo runs, predict that M, S and 
C stars should occur in the ratio .65:. 04:. 31. This is a 
(TP-AGB) C/M star ratio of 0.47, in rough agreement 
with observations - provided the observed ratio is indeed 
about 0.6; the LMC C/M star ratio on the TP-AGB is 
a rather uncertain quantity. 

The observed LF of C stars in the LMC was repro- 
duced from Costa & Frogel (1996). Fig. |^ gives exam- 
ples of synthetic LFs that approximate the observed LF. 
The synthetic LFs were smoothed with a cr = 0.2 gaus- 
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sian, as in GJ93. Structure on small scales would not be 
picked up in observations anyway; there is considerable 
scatter about the bolometric correction law (Costa & 
Frogel 1996, their Fig. 1). The distance modulus to the 
LMC was taken as 18.5. 

As a note of caution, Marigo et al. (1999) have 
shown how the C star LF is more sensitive on metal- 
licity than it is on either the star formation history or 
the IMF. The metallicy dependence of the prescriptions 
we use also results in such a sensitivity, but an extrap- 
olation from the LMC to the solar neighbourhood may 
seem far-fetched. It is assumed that none of the free 
parameters are metallicity dependent. 
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A2.2 The initial-final mass relation 

The initial-final mass for the 'standard' galactic model 
is shown in Fig. A2, where it is compared to the data 
from Weidemann & Koster 1983 and Weidemann 1987 
(data points from FC97). Final masses are derived from 
the surface gravity, or the stellar radius; there are thus 
sometimes two data points per object, with the same 
Mi. The agreement is satisfactory, except at low ini- 
tial masses, where the model initial-final mass relation 
might be an overestimate (it should be kept in mind 
that the 'observed' initial-final mass relation is in fact 
model dependent). Wagenhuber & Groenewegen (1999) 
reach the conclusion that the GJ93 approximations for 
the core masses at the first thermal pulse (which are 
used here, and largely determine the final mass), over- 
estimated most other model results at low initial mass, 
but the discrepancy is small. The scatter in the model 
Mi — Mf stems from the metallicity dependence of the 
analytical prescription for the core mass at the first 
thermal pulse (it is smaller for the prescription in Wa- 
genhuber & Groenewegen 1999). 

A3 Evolution on the TP-AGB for 
characteristic initial masses 

The synthetic AGB evolution of models with {Mi, Zi) = 
{(1.5, 0.01); (5, 0.02)} are shown in figure ^ The case 
with Mi — 1.5 Mq shows how towards the low mass 
end of C-rich PNe progenitors, a low number of thermal 
pulses gradually builds up a ^^C surface overabundance. 
With Mi = 5 M0, the first ~ 2 10^ yr correspond to sur- 
face C enrichment, until core growth is high enough for 
Tb to reach CNO processing values. The total number 
of thermal pulses reaches 143, and at i = 10^ yr the en- 
velope mass drops below the adopted critical mass for 
the quenching of HBB. 



Figure A2: The initial- final mass relation from the 
galactic model in this article, and as observed. Dia- 
monds and open-circles correspond to final mass esti- 
mates based on surface-gravity or stellar radius, respec- 
tively. 
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Figure A3: Synthetic evolution for characteristic initial 
masses. The results for the surface ^^C/^'^C, log(N/0), 
C/0, the mass loss rate and bolometric absolute mag- 
nitude are shown from top to bottom, for initial masses 
1.5 and 5 M0 from left to right). 
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